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Abstract 

In this paper changes in wind speed and wind direction from a measured wind field 
are being analyzed at high frequencies. This is used to estimate changes in the angle 
of attack (AOA) on a blade segment over short time periods for different estimated 
turbine concepts. Here a statistical approach is chosen to grasp the characteristics of 
the probability distributions to give an over all view of the magnitude and rate of the 
changes. The main interest is the generation of basic distributions for the calculation 
of dynamic stall effects and stall flutter due to wind fluctuations. 
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1 Introduction 

Dynamic stall is a concept to calculate strong, yet hard to predict, load fluctuations on 
wind turbines. One of the main causes for dynamic stall and stall flutter on wind turbines 
are sudden changes in wind speed and wind direction leading to changes in the angle of 
attack (AOA) on wind turbine blades. Both angle of attack and dynamic stall are part of the 
2D calculation concept used by most blade element moment (BEM) or lifting line model 
theory programs to calculate the main aerodynamic properties for wind turbines HlGl- 
Dynamic stall is mostly taken into account by empirical models based on the knowledge of 
the amplitude and fequency of the changes in the angle of attack [3]. 
Most of the research on the dynamic stall models evolved from helicopter research like the 
ONERA model [i9j or the Beddoes-Leishman model ifTOll and were later transfered to wind 
turbine applications. At helicopters dynamic stall mainly appears due to a rapid change in 
AOA induced by a - from the rotor perspective - yawed inflow IfTTI caused by the flight 
velocity. The changes in AOA and its rates of change are periodic in time and thus can 
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easily be calculated. Such change rates are also often used for wind turbine dynamic stall 

modeling ma ini. 

In case of a natural wind field the situation is different, as changes in the AOA are induced 
by the incident wind field which appears in most cases to be highly turbulent. To calculate 
the amplitudes and frequencies of changes in AOA the distributions of changes in wind 
direction and wind speed have to be investigated. 

The fluctuations of the wind has attracted attention in research in recent years (e.g. Il4l ll5lll6l ) 
The investigation on the influence of such fluctuations on the AOA is a relevant conse- 
quence. Bierbooms and Veldkamp did some research on wind speed fluctuations and the 
resulting loads on wind turbines |[7l|[8||. 

For the load estimation on wind turbines, wind field generators using a wind field based on a 
Gaussian distributed wind speed fluctuations are applied quite often (e.g. IITTll IfTSlI [|T9lll20ll ). 
This Gaussian approach leads to an underestimation of extreme values of the fluctuations 
in wind speeds (see e.g. fTDi ). 

Changes in wind direction lead to a yawed inflow condition, which is one reason for wind 
turbines to face dynamic stall effects [[T4lllfT5ll .The analysis of wind directional changes has 
been described by van Doorn et al. lfT6ll . The methods are picked up in this paper and 
applied to a measured wind field for further statistical analysis. 

If the fluctuations in wind speed and changes in wind direction would be strongly corre- 
lated, it would be straight forward to develop a comprehensive analytical model for the 
fluctuations in the AOA. However this does not often seem to be the case [|22ll . Therefore 
we propose a direct calculation of the fluctuations in the AOA from changes in both, the 
wind speed and the wind direction. 

Due to the unsteadiness of the wind field and the response of the wind turbine, the actual 
AOA is hard to determine. Different approaches to resolve this problem have been pro- 
posed ll23l ll24l ll25l ll26l . In this paper we assume that the actual AOA is unknown. Yet the 
changes in the AOA due to changes in wind speed and wind direction are evaluated. We do 
so by estimating the changes in the AOA over short period of time. Therefore will assume 
certain simplified models for such changes on an arbitrary chosen position on a rotating 
blade. 

Here we will explain the methodology to calculate the local changes in the AOA, using a 
measured wind field as input. At first a short analysis of the wind field data will be done. 
Next, for two different wind turbine modeling principles, the changes in AOA are calcu- 
lated and analyzed. Included is an analysis of the changes in AOA only due to changes in 
wind direction. In the end an analysis of the time scales on which the changes appear and 
the characteristics of their distribution are presented with the aim to improve the investiga- 
tion of dynamic stall effects and their modeling for wind turbines. 



2 Characteristics of the measured wind 

Before we present a method to calculate the changes in AOA, the measured wind field we 
used for the demonstration shall be discussed here on the bases of its statistical aspects. 
The data from a measurement campaign at a wind farm at Meerhof near Paderborn in Ger- 
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many is being used throughout this paper [|27ll . The data was measured using a Gill R3-50 
ultra sonic anemometer in 98m height with a sampling rate of 50 Hz. If not otherwise 
stated, it has been averaged to a 10 Hz sampling rate. Using the Taylors frozen turbulence 
hypothesis this sampling rate corresponds to a length scale of 1 m at a wind speed of 10 
m/s. Wind coming from the sector from 40.5° to 133.5° has been excluded from the anal- 
ysis, since a wind farm was near in these wind directions. This way an influence by the 
wake of the wind farm on the wind measurement could be avoided. For the investigation a 
dataset of seven consecutive days from December 1st to 7th, 2005 has been used, leading 
to a dataset of 6 ■ 10^ samples at 10 Hz. 

To give an impression of the wind field used in this study we describe different relevant 
characteristics of the wind field and the methods for the characterization of turbulent fluc- 
tuations. 

The mean measured wind speed over all data was w = 6.96 m/s with a turbulence intensity 
of / = 13.3% according to IEC-6 1400-1 Standard 11281 . calculated from the average stan- 
dard deviation of 10 minutes periods. 

For the further calculations some approximations have been made. We exclude wind speeds 
of |m| < 2 m/s in the analysis. Changes in wind direction at low wind speeds are very sen- 
sitive to measurement errors. Very low wind speeds are also irrelevant for wind turbines, 
due to the lack of induced relevant loads. 

Since the dataset consisted only of one point measurements, we will also make the assump- 
tion for the calculations that the wind field and its turbulence is homogeneous. Even if this 
is in reality not the case, implementing a different wind field on the method will be straight 
forward. 

For the research of dynamic stall we are mainly interested in the extend of the changes at 
certain periods. Therefore the wind field and its turbulent structure was further analyzed 
using incremental statistics. They are a good measure for fluctuations at steps x. Increments 
of a property ^,{x,t) are defined as 

^^{x)=^{x,t + x)-^{xj) . (1) 

where x is the time scale to be considered. The statistics of the increments allow a time 
scale resolved analysis of the fluctuations of a magnitude ^. 

For the measured data the probability density function (pdf) of the wind speed increments 
for a time scale of x=10 Hz is given in fig. [T]in semi-logarithmic presentation. We observe 
a typical intermittent behavior of the pdfs, which has been described e.g. by Bottcher |[2T]| . 
The distribution deviates strongly from a fitted Gaussian curve. In our case the empirical 
probability of wind speed changes at 10 Hz and of magnitude 0.85 m/s is already by a 
factor > 10"^ higher than a corresponding Gaussian distribution would predict. For larger 
increments this deviance increases by further orders of magnitude. 

The deviance from a Gaussian curve can be quantified by the kurtosis and the skewness, 
where the kurtosis can be used to describe a symmetric deviation from a Gaussian curve 
and the skewness to describe the asymmetry of the curve. In this paper the excess kurtosis 
is taken into account given by 
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Figure 1: Probability density function (pdf) of all wind speed increments for ^ = 10 Hz. 
The typical intermittent behavior is seen by the strong deviations of the tails from a fitted 
Gaussian distribution. Here like for all other plots we used only wind speeds at \u\ > 2 m/s. 



The excess kurtosis is normed so that y = for a Gaussian distribution and y > for inter- 
mittent pdfs. 

The skewness is calculated by the standardized third moment: 

V = ^^'=i(^'~^)^ (3) 
It follows that V = for symmetrical pdfs. 

For a characterization of the wind field, we can determine the pdfs of the wind speed 
increments at different time scales x (see fig. H]). In fig. [2]b) the corresponding values of the 
excess kurtoses are shown as a function of time scales between 0.1s and 3s is shown. With 
an increase in time period of wind speed increments y becomes smaller, which means that 
the shape of the pdfs approach slowly the shape of a Gaussian curve. A linear regression 
indicates a characteristic y ^ x^'^-^^ An extrapolation would lead to a kurtosis of a Gaussian 
curve at a time scale of 150 s. This is however a time scale far above the scales at which 
most control systems adapt wind turbines to the surrounding conditions. Thus it is of no 
further interest for our considerations here. 

So far we have look just at the wind speed \u\ and its increments. For our purpose it is 
also important to take the wind direction (|) and its increments into account. 
The pdfs of the increments of the wind direction show again an intermittent form (see fig. 
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Figure 2: pdfs wind speed increments for different time scales x on a semi-logarithmic 
scale a). The values of the excess kurtosis of the pdfs give an idea of the magnitude of the 
deviation from Gaussian distributed pdfs with a y = 0. 




Figure 3: a) probability density functions of wind direction increments at different time 
scales on a semi-logarithmic scale and b) the excess kurtoses. 
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[3]). The values of the excess kurtoses are lower than for the wind speeds. Fig. [3] makes 
clear, that wind direction changes are large and appear rapidly. Changes of 10° and more 
within 0.1 seconds do occur. For larger time scales, still less 1 second, changes of up to 
20° occur. The double logarithmic presentation of the excess kurtoses in fig. [3]b) leads to 
a similar power law with y x"*^-^^ for scales 0.5 < x < 3.0. An interpolation leads to a 
kurtosis of a Gaussian curve a.tx^A.6s. 

As both changes in wind speed and in wind direction can cause changes of the angle 
of attack, the correlation of the appearances of both phenomena is of particular interest. A 
possible measure to quantify this is the so called gust directional index (GDI) 

^ abs{u{t + x)-u{t)) ^ abs{(^{t + x)-(^{t)) 
max{abs{u{t + x)~u{t))) max{^{Dir{t + x) - ^{t)))' 

as described by Hansen et al. If22]| . Here u{t) is the wind speed at time t, x is the time scale 
and (|)(?) the wind directione. The GDI E [0,2] takes into account two simultaneous incre- 
ments and <^i{t), with GDI{x) = 2 if both magnitudes show extreme events occurring 
in a correlated way and GDI{x) = for complete non-correlation. For quantification the 
probability of strongly correlated events is regarded [|22]| . 

For an 1 Hz sampling rate the probability of the GDI > 1 .98 for x = 10 s was at 1 .3 x 10—5. 
For a sampling rate of 10 Hz this value for x = 10 s decreases to 1.5 x 10^^. It seems a 
large simultaneous change in wind speed and direction is not likely to appear, the numbers 
of large GDI events are comparatively low. 

The low values indicate for the question of the changes in angle of attack, that changes in 
wind speed and direction could mostly be considered separately for this wind field. Never- 
theless this does not seem to be the case for all wind situations (compare ll22ll '). 



3 Analysis of fluctuations in the AOA 
3.1 The calculation concept 

Instead of focussing on aspects of the wind field components separately we proceed to 
calculate the changes in AOA by a comprehensive approach. To do so a model of an 
estimated wind turbine is need. The real AOA on wind turbines is strongly dependent on 
its control system. Here we do not intend to go into details of a specific control system. 
Instead we consider two different simplified cases: 

a) A classical fixed tip speed ratio turbine that adapts tip speed and pitch angle to the 
changes in the averaged wind of a certain period x^ by a fast controlling system. We 
will take here x^ = 2 s. 

b) A turbine with a fixed rotational speed, that adapts the pitch angle to the averaged 
wind of the last 2 seconds. 
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The models will further be labeled as case a) or b) respectively. 

For all cases of the models we find a finite time Xa for the control system to pose a wind 
turbine in its optimal working condition. In case of wind changes on time scales x < Xq 
the control system is not able to react quickly enough. For these situations the fluctuations 
of the wind conditions e.g. grasped by increments are of importance. This aspect will be 
studied in the following. 

The used estimations of turbines enable us to determine a change in AOA by a change in 




rpm 



Figure 4: The changing AOA can be calculated by the law of cosines. In grey is the original 
wind velocity and direction - black is the new one after a time x. 

wind situation. Such changes can then be calculated by using the geometry shown in fig. 
|4l where 



A4) = 4)(/'+x)-4)(0. 



(5) 



Here <^{t) is the averaged angle between the incoming wind and the rotor plane determined 
by the "slow" yaw system. 

The general equations for the changing AOA Aa in fluctuating winds are, in case a) of the 
fixed tip speed ratio turbine 



Aa = arccos 



u{t)^-lu{t + x) sin(A(^) 



'{u{t)^)^ + {ju{t + x))2 - ^u{t + x)u{t) ^ sin(A^) 
For the fixed rotational speed turbine case b) we use 



'■rpm 



60 



(7) 



where n,- is the rotational speed in number of rounds per minute. For turbines with a rotor 
radius of i? = 40 m, typical values of Ur lay between 10 and 20. Thus we obtain: 



Aa = arccos 



Urpm -\u{t + x) sin(A(|)) 



^"rpm + + T^))^ + T^)«rpm sin(A(|)) 



arccos 



'■rpm 



'u}pm + {l<t)Y 



(8) 
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Here in both cases u{t) is the absolute value of the wind velocity at the time t, is the tip 
speed ratio, which was taken to be 7 for case a). R is the total rotor radius, which is con- 
sidered to be = 40 m in this paper, r is the position of the blade, which is to be regarded. 

'-) 

The factor of | in both equations reflects the estimated deceleration of the wind due to the 
blockage effect. 

All calculations were done at a position of half of the blade length r = 20 m. Although the 
absolute changes in AOA depend on r, the form of the distributions of the increments do 
not. 



3.2 Overall estimation of changes in AOA 

In fig. [Hand fig. [3] we can conclude that the timescales of rapid fluctuations do not differ 
very much between wind speed change and wind direction change. Further we can con- 
clude from the low value of the GDI in section [2] that we cannot expect a strong correlation 
of the fluctuations in wind speed and wind direction. Therefore it makes sense to investi- 
gate the changes in the AOA by regarding the over all change in wind speed and direction 
at the same time, calculated by equations Q and ([8]). 

Here the models for case a) (const, tip speed) and case b) (constant are taken into 
account. We use la = 2 s. Fig. [5]presents the pdfs of the increments Aa for both cases, in 
a) with 

Uypm — 20 in b) Urpm — 10 representing the limits of the range of rotational speed 
for a typical R=40 m turbine. To emphasize the difference in shape of the pdfs between the 
two cases a) and b) the results have been plotted together for each M,p;7i- 
The pdfs in fig. [5] for both assumed models of the turbines show a strong intermittent 
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Figure 5: Pdfs of the changes in AOA due to change in wind over the period of = 2s. In 
both figures a) and b) the changes of the calculations for constant rotational speed against 
(case b)) constant tip speed ratio (case a)) are given. In figure a) the const, rotational speed 
is given by a rotation of Urpm = 20, whereas in fig. b) rirpm = 10. Note that the pdf of case 
a) with depicted by the round symbols is the same in both plots. 
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Type 


Skewness 


Kurtosis 


a 


= 7 


0.20 


1.8 


0.59° 


^rpm 10 


0.16 


3.1 


0.71° 


^rpm 20 


0.17 


3.2 


0.36° 



Table 1: Characterization of the pdfs for the change in AOA due to wind fluctuations - 
including wind speed and direction. 



structure. In case of the case a) turbine with constant tip speed ratio, the intermittency is 
less distinct. This is also quantified by the kurtoses of the pdfs in table [TJ The kurtosis of 
the case a) pdf is much lower yet the standard deviation a is relatively high compared to 
the pdfs of case b). The reason lies in the adaption of the case a) turbine rotation to the 
wind speed filtering some of the extreme fluctuations. 

Further a slight asymmetry in the pdfs could be observed which are quantified by the skew- 
ness in table [U Since the fluctuations in wind speed should not cause such an asymmetry, 
this effect should be a result from the changes in wind direction. To quantify the influence 
of the changes in wind direction on the pdfs of the changes in angle of attack, we will pro- 
ceed to analyze this effect separately. 



3.3 Maximum changes of attack due to sudden yaw angle changes 

In the following we analyze the effect of the changes in wind direction on the changes in 
angle of attack. 

We regard the maximum changes in AOA under the following conditions: The wind speed 
is expected to be constant, the averaging period x^, = 2 s and = 1 for case a) and rir — 
10, 15, 20 rpm for the case b) model. Thus for the = 15 case a rotor blade would perform 
half a revolution. 

Using these assumptions the model equations ^ and ([8]) change to 




^-|sin(A(^) , , - - . 

Aa = arccos | k i ^ ^ | _ Cretan | -— | , (9) 

;^)2 + (f)2-f^sin(A^), 
for turbines with constant tip speed ratio and to 







1 — arctan 1 - 






Irr j 




R / 



Aa = arccos I — rpm 3 w v 1 _ ^^^^^ j — I (10) 



^rpm + (f "(0)' - lu{t)u,p,„ sin(A(^) / y Jujp,, + (f 



for a turbine at a constant rotational speed. 

Using A(|)(Xa) we obtain the pdfs for the changes in AOA presented in fig. [6l 

The pdfs for the AOA in fig. [6] show a remarkable difference between the case of a 
turbine running at constant tip speed ratio and the turbine running at constant rotational 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 

da [°] da [°] 

Figure 6: Pdf of the changes in AOA due to change in wind direction within 2 s. a) depicts 
the case a) of the turbine with a constant tip speed ration. The figures b), c) and d) show 
the pdfs for the case b) model rotating at a constant speed of rirpm = 10, 15 and 20 rpm 
respectively. 
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Type 


Skewness 


Kurtosis 


a 


= 7 


-1.7 


7.5 


0.28° 


^rpm 10 


-0.9 


8.2 


0.23° 


^rpm =15 


-1.6 


10.9 


0.11° 


^rpm — 20 


-2.2 


14.1 


0.06° 



Table 2: Characteristics of the pdfs of Aa due to change in wind direction 

speeds. Due to the constant ratio between the wind speed and the rotational speed the pdf 
of the case a) model shows a distribution limited towards positive deviations. While the 
pdfs for the case b) model show an intermittent curve, which is strongly asymmetric. The 
base functions ^ and (flOl) lead to a higher scatter toward negative angles in all cases. As a 
result all pdfs in the graphs have a negative skewness and very high values for the kurtoses 
given in table [2l Comparing the characteristics of the pdfs in table [2] to the values of the 
pdfs obtained for the complete wind field in table [B we conclude that the influence of the 
changes in wind direction on the changes in AOA are small. 



4 Deeper analysis of the statistics of Aa 

To develop models for the fluctuations in the AOA, the rate of the changes over time have to 
be regarded. This helps to derive a model for the turbulence of the changes in the AOA and 
as a result realistic estimations of extreme values. Also the rates of Aa over time is of major 
interest for dynamic stall modelling. Therefore the fluctuation of an estimated AOA at dif- 
ferent time scales were calculated. To gain a quantative image of the magnitude of the time 
scales of fluctuations in angle of attack, instantaneous changes in the AOA without averag- 
ing have to be regarded. Therefore here a modified (case c)) version of equation ^ is being 
used for the further calculations. Thus only a comparison between the current wind and the 
one at a time different x is done, so that and u{t) have to be replaced by (|)(?) and u{t) 
in equations © and ([8]) respectively. The time scale was henceforth varied from 10 to ^ Hz. 

It can be seen in fig. |7] a) a main portion of the changes in Aa appear to evolve at short 
time periods. For the characterization of the distributions we plotted the kurtosis over the 
time scales in fig. |7]c). It can be observed that the kurtosis of the pdfs shrinks only slowly 
for time scales of 10 - ^ Hz. As an example of the deviation from a Gaussian distribution 
fig- |7]b) depicts the distribution for a time scale of x = ^ s. 

The knowledge of the pdfs make reconstructions of such distributions by analytical models 
in case of wind speed distributions possible ll29ll . If such models are applicable also to the 
change in AOA Aa, they can further be used to optimize models for dynamic stall on wind 
turbines. Castaing et. al. proposed to fit the pdfs of Am to a function that would for Aa look 
like 

N 1 / Aa2\ / ln^(o/oo)\do 
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Figure 7: The pdfs for different time scales are given in a). As expected Aa tends to be 
broader for longer time scales. The pdfs remain nevertheless intermittent all the time. Fig. 
b) depicts the pdf of the jq s time scale against a Gaussian curve. Both figures a) and b) 
are calculated for a fixed rotation of nrp,„ = 20 rpm. The rate of the intermittency over 
the time scale is given by the kurtosis in fig. c) for fixed rotations of rirpm = 10 rpm and 
firpm = 20 rpm respectively. Fig. d) shows the broadening of pdfs by its standard deviation 
and the time scale x. The line in d) depicts the slope of a linear function of 4 for the case 
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where Gq is the most probable variance of Aa OOll . By determining the X we thus gain the 
analytical fit function P{Aa^) of the pdf for the time scale x. Since X'^ can be taken from 

X2oc/„(Y), (12) 

the knowledge of the different statistical properties as o(x) and y(x) is essential, they de- 
scribe the pdfs at different conditions. In the following they shall be described in more 
detail to give an overview over the characteristics of the pdfs. 

The distribution of the pdfs given in fig. |7]a) only broaden at a mediocre rate as the time 
scale increases. In fig. |7]d) the double logarithmic plot the standard deviations a of the pdfs 
of Aa over the different time scales are shown. For time scales x > 0.6 s the broadening of 
the standard deviation can be well described by a power law of a ~ X3. If we assume the 
Taylor hypothesis to be valid in that range 

Ar =<u> At, (13) 

the slope of the curve in fig. |7]d) coincides with the turbulence theory of Kolomogorov 
from 1941 for isotropic and homogeneous turbulence OTI in case of the second structure 
function, stating 

S:ir)ocr'i. (14) 

Here S is the structure function of the n-th order. As the variance is the structure function 
of the 2-nd order, o °^ could be expected for Gaussian turbulence. So even though the 
overall pdfs for the change in angle of attack for time scales 0.6 < x < 3.0 s show a non 
Gaussian behavior, the relations for the structure function for these time periods behave 
according to ll3T1l like a Gaussian distribution. 

A possible reason for this could be, that the overall distribution consists of a superposition 
of many mainly almost Gaussian distributions. In fig. [8]the pdfs and the kurtoses over time 
scales are shown for wind speeds under the restricting condition 6 < |m| < 8 m/s. Already 
the pdfs can clearly be identified to be a lot closer to a Gaussian shape then it was the case 
for the non-conditioned wind speed in fig. |7^). Though x < 2 s the kurtosis of the pdfs is 
still positive, it is a lot smaller than for the over all wind field. This explains why obviously 
equation (fT4l) is in our case still valid. 

As could be seen in fig. |71 for time scales x < 0.6s the relation of equation (fT4l) does not 
hold anymore. This is also reflected for the conditioned pdfs in figl8]for small time scales 
X < 0.5 s. The kurtoses do not decrease or approach a value close to y = 0. 
Using this background it is a reasonable attempt to fit a reconstructed pdf to the data to gain 
an analytical function with equation (fTTI) for the distribution. Such a fit has been done for 
the conditioned pdfs in fig. [8] a). The pdfs give to a good extend the characteristics of the 
changes in the angle of attack at a certain position of a wind turbine blade within a given 
period of time. 

Since for the calculation of dynamic stall not only the magnitude, but also the time in 
which changes in angle appear are important, the rates of the changes in time are being 
presented here. It is straight forward to calculate the distribution of da/dt for a sampling 
rate of x = 10 Hz. Fig. [9] a) shows the distribution for Aa/x for a rotation of nrpm = 20 
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log(T) [s] 



Figure 8: Looking at pdfs of Aa under certain wind speed conditions. In a) pdfs for differ- 
ent time scales under the condition of 6 < \u\ < 8 m/s show increasingly Gaussian shapes 
the larger the time scale gets. This is quantified by the plot of the kurtoses in b) for different 
conditions. Where the original curve shows way larger values. 

rpm. In fig. [9] the rates are plotted for the standard deviation of the pdfs and the maximum 
values for different x. 
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o const rpm = 20 
o Gauss Kurve 



b) 
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Figure 9: Probability density function of Aa/x for x = 10 Hz graphed with a Gaussian 
distribution for comparison, giving an idea for da/dt a). The maximum rates and the rates 
for the standard deviation of the distribution for different x is given in b). 

5 Discussion and Conclusion 

Inspired by the strong simplifications in the modeling of dynamic stall, we have shown an 
analysis of the statistical properties of the changes in the angle of attack estimated from 
the wind field. After a short analysis to characterize the wind field, we proposed a simple 
comprehensive model to estimate the changes in the angle of attack on a turbine blade at a 
certain rate of time. Basic turbulence characteristics from the wind field were thus trans- 
formed into characteristics which could be used for the modeling of dynamic stall on wind 
turbine blades. 

The statistical properties of the wind speed and the wind direction showed the typical in- 
termittent behavior. Even though the latter had a tendency to be less intermittent. Both 
changes in wind cause changes in the angle of attack on blades of a rotating turbine. Such 
changes therefore reflect the distributions of the changes in wind and showed in general 
also an intermittent distribution. Different cases were regarded in the study. 
Two model cases were used to calculate the changes in AOA: case a) for constant tip speed 
ratio and case b) constant rotational speed. For both cases the changes in AOA with a mov- 
ing average at a time scale x^ = 2 s were analyzed. In case a) the adaption of the rotational 
speed to a moving average of the wind speed caused the pdfs of the changes in AOA to be 
less intermittent than the ones for turbines running at a constant rotation rate. 
To analyze the asymmetry in the pdfs, the same cases have been evaluated for the instance 
of a change in wind direction only. Here pdfs for the changes in AOA show a strongly 
skewed behavior. However for the over all pdfs of the changes in AOA caused by wind 
speed and wind direction changes, it could be seen that the changes in wind speed domi- 
nate the effects on the changes on the AOA. 

Since the aim was to find models for the changes in AOA that could be used for dynamic 
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stall modeling the turbulent structure of the changes have been further analyzed. It was 
found for turbines running at a constant speed, that for time scales x > 0.6 s the standard 
deviation of the pdfs broadens by a power law ^ x^. From this was concluded that the 
intermittent distributions over the complete time consist most likely of a superposition of 
mostly Gaussian distributions. For time scales x < 0.6 s this is not true anymore. So for 
AOA changes at time scales of x < 0.5 s larger values of the kurtoses and smaller values 
of the standard deviation are to be expected - indicating intermittent, non-Gaussian fields. 
Due to the small time scales, this is expected to lead to local aerodynamic effects on the 
blade as dynamic stall and will not happen on scales to which an existing controlling sys- 
tem would be able to react. 

Since the main interest was to find characteristic distributions for the dynamic stall mod- 
eling the resulting conditioned pdfs for the changes in AOA have been fitted. The results 
showed slight deviations caused by the changes in wind direction. However the major char- 
acteristics were grasp very well. This makes new approaches for the modelling of dynamic 
stall for the optimization of wind turbine design possible. 
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